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ABSTRACT 

We demonstrate the existence of a secular back-reaction on inflation using 
a simple scalar model. The model consists of a massless, minimally coupled 
scalar with a quartic self-interaction which is a spectator to A-driven in- 
flation. To avoid problems with coincident propagators, and to make the 
scalars interact more like gravitons, we impose a covariant normal ordering 
prescription which has the effect of removing tadpole graphs. This version of 
the theory exhibits a secular slowing at three loop order due to interactions 
between virtual infrared scalars which are ripped apart by the inflating back- 
ground. The effect is quantifled using an invariant observable and all orders 
bounds are given. We also argue that, although stochastic effects can have 
either sign, the slowing mechanism is superimposed upon them. 
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1 Introduction 



The application of a force field in quantum field theory generally rearranges 
virtual quanta and thereby induces currents and/or stresses which modify 
the original force field. This is the phenomenon of back-reaction. Famous 
examples include the response of QED to a homogeneous electric field [Q 
and the response of generic matter theories to the gravitational field of a 
black hole 0. In the former case, virtual e~^e~ pairs can acquire the energy 
needed to become real by tunneling up and down the field lines. The newly 
created pairs are also accelerated in the electric field, which gives a current 
that reduces the original electric field. The event horizon of a black hole 
also causes particle creation when one member of a virtual pair passes out of 
causal contact with the other by entering the event horizon. As the resultant 
Hawking radiation carries away the black hole's mass, the surface gravity 
rises. 

Parker was the first to realize that the expansion of spacetime can 
lead to the production of massless, minimally coupled scalars. Grishchuk 
later showed that the same mechanism applies to gravitons. Production of 
these particles is especially efficient during inflation because virtual infrared 
pairs become trapped in the Hubble flow and are ripped apart from one 
another. Since the mechanism requires both effective masslessness on the 
scale of inflation and the absence of conformal invariance, it is limited to 
gravitons and to light, minimally coupled scalars. 

Our special interest is the back-reaction from this process. We believe 
that the gravitational attraction between virtual infrared gravitons gradually 
builds up a restoring force that impedes further inflation. This mechanism 
offers the dazzling prospect of simultaneously resolving the (old) problem 
of the cosmological constant and providing a natural model of inflation in 
which there is no scalar inflaton . The idea is that the actual cosmological 
constant is not small and that this is what caused inflation during the early 
universe. Back-reaction plays the crucial role of ending inflation. 

The purpose of this paper is to establish that there is a significant back- 
reaction. This might appear obvious in view of the limitless extent of particle 
creation otherwise. It is easy to show that about one infrared pair emerges 
per Hubble time in each Hubble volume. Denying that there is significant 
back-reaction implies that an observer can watch this go on in the space 
around him for an arbitrarily long period without feeling any effect. 
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However, the preponderance of expert opinion is highly doubtful about 
the existence of a significant back-reaction. Some are concerned with the 
methodology of previous studies of back-reaction. In the pure gravity model 
described above, the expectation value of the gauge-fixed metric was com- 
puted in the presence of a state which is free graviton vacuum at t = 0, and 
the resulting invariant element was reported in co-moving coordinates, 

{fl \g^^{t, x)dx^'dx''\ Q) = ~dt^ + e^^^^dx ■ dx . (1) 

When two loop effects are included the expansion rate is. 
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+ 0(G^) , (2) 



where G is Newton's constant, A is the cosmological constant and H = 
is the Hubble constant of the locally de Sitter background 0. Unruh has 
criticized the procedure of taking the expectation value of the metric first and 
then forming it into an invariant measure of the expansion rate Linde has 
no objection to gauge fixing but he believes that expectation values obscure 
important stochastic effects 

There are also concerns about the putative physical mechanism behind 
screening. Some doubt the causality of gravitational interactions between 
particles which have been pulled out of causal contact with one another. 
Others concede the reality of such a residual interaction but maintain that it 
must be redshifted into insignificance by the inflationary expansion of space- 
time. Finally, there are those who insist that back-reaction must degenerate 
to a small increase in the expansion rate because the time and space aver- 
age stress-energy tensor induced by inflationary particle creation is that of a 
small, positive cosmological constant. 

Our response to plausible methodological concerns is cooptation. We 
believe that a physical process such a back-reaction will manifest itself in 
any reasonable formalism. To address Unruh's objection we form the metric 
operator into an invariant measure of the local expansion rate before taking 
its expectation value To address Linde's objection we set the formalism 
up so that stochastic samples of this expansion operator can be taken instead 
of expectation values [|10] . 



These changes are too complicated to be implement in the pure gravity 
model but they can be carried out in a scalar analog which possesses the 



2 



same combination of attractive self-interactions between massless, confor- 
mally noninvariant quanta. The model is essentially a massless, minimally 
coupled scalar with quartic self-interaction, 

+ counterterms + ordering corrections , (3) 

in a locally de Sitter background. To make this system more similar to 
gravity we remove tadpole graphs through a procedure of covariant normal 
ordering that preserves conservation of the stress-energy tensor. When this 
is done one can apply the same formalism which was used for pure gravity 
to show that the expansion rate is slowed by an amount which eventually 



becomes non-perturbatively strong 

h{t)=H\^l-^^[{HtY + 0{HH')\+0{\\G^)^ . (4) 

We shall demonstrate that this result is not changed by using an invariant 
operator ala Unruh, and that it is superimposed upon an indeterminate effect 
of order A when one takes stochastic samples ala Linde. 

Section 2 addresses the various non-methodological objections (causality, 
redshift, and modeling the source as a homogeneous, classical fluid). Section 
3 motivates our decision to order the scalar analog theory so as to subtract 
off tadpoles. Section 4 explains the actual procedure for accomplishing this. 
In Section 5 we show that the expectation value of the invariant expansion 
operator agrees exactly with (H). We also argue that, while a stochastic 
sample can show an effect of either sign at order A, there is no significant 
change in the order A^ result. Our conclusions comprise Section 6. 



2 Physics of inflationary back-reaction 

The purpose of this section is to review the physics of inflationary back- 
reaction so as to answer the three non-methodological objections which were 
summarized before. We shall also give a simple explanation of why back- 
reaction induces an ever-increasing, negative vacuum energy in perturbation 
theory. The aim here is not rigor — that is supplied by the detailed cal- 
culations of Section 5. We seek rather to explain in simple terms why the 
calculations turn out as they do. 
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It is worthwhile recalhng that the inflationary production of gravitons 
(and light, minimally coupled scalars) is a straightforward consequence of 
the Uncertainty Principle, the existence of a causal horizon during inflation, 
and the simultaneous masslessness and absence of conformal invariance of 
the quanta being produced. The Uncertainty Principle requires all quantum 
degrees of freedom to possess 0-point motion. In a quantum field theory this 
means that virtual quanta continuously emerge from the vacuum. Although 
the fact of their emergence does not depend upon the background geometry, 
what happens to them subsequently does. In particular, if there is a causal 
horizon then long wave length quanta emerge out of contact with one an- 
other and so can never recombine. For this to happen with an appreciable 
amplitude the quanta must of course be massless on the scale of the horizon. 
They must also be sensitive to the local geometry. For the conformally fiat 
backgrounds characteristic of long periods of inflation this implies that the 
quanta must not possess classical conformal invariance. That is why gravi- 
tons and light, minimally coupled scalars experience super-adiabatic amplifi- 
cation during inflation whereas vector gauge bosons, spin 1/2 fermions, and 
conformally coupled scalars do not. 

In considering the gravitational back-reaction from inflationary particle 
production note first that there is no buildup of particle density because 
the 3-volume expands as new particles are created so as to keep the density 
constant. When a new pair is pulled out of the vacuum the one before it is, on 
average, already in another Hubble volume. However, the gravitational field 
is another thing. The created particles are highly infrared (by the standards 
prevailing at the time) so they do not carry very much stress-energy, but they 
do carry some. This must engender a gravitational field in the region between 
them. Because gravity is attractive, this field must act to resist the Hubble 
fiow. This is a very small effect because gravity is a weak interaction, even 
on the scales usually proposed for inflation. The feature that permits it to 
become significant is the cumulative nature of the effect. Even after a newly 
created pair has been pulled into distinct Hubble volumes its gravitational 
field must remain behind to add with those of subsequent pairs. If nothing 
else supervenes to end infiation first the gravitational field must eventually 
become nonperturbatively strong. 

There is little doubt that infiationary particle production goes on be- 
cause it is the usual explanation for the primordial spectrum of cosmological 
perturbations [O, IT3| whose imprint on the cosmic microwave background 
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has been so clearly imaged by the latest balloon experiments Nor is 

there any real doubt that some gravitational back-reaction accrues from the 
process, since the usual theory of structure formation holds that we live in 
complex structures resulting from the gravitational collapse of primordial 
fluctuations over the course of ten billion years. The real issues are, whether 
or not there is a gravitational response during inflation, whether or not this 
response grows with time, and whether or not the response slows inflation. 

Those who argue against any response at all base themselves on causal- 
ity. The source for our process is virtual particles which rapidly fall out of 
causal contact with one another. Therefore, how can they interact in any 
way? The flaw in this argument is the picture it implies of exchange forces 
being maintained by instantaneous interaction across some arbitrary surface 
of simultaneity. Were this correct an outside observer could not feel the grav- 
itational attraction from a pebble after it had entered the event horizon (in 
some arbitrary frame) of a black hole. In fact exchange forces are maintained 
by virtual quanta which propagate causally, and those which carry the grav- 
itational force from an object falling into a black hole originate in the region 
outside the event horizon. The same is true of objects which fall out of causal 
contact in an inflating universe: they continue to feel a gravitational force 
from one another that is carried in a completely causal manner from far back 
along their past light cones. The process is depicted in Fig. 1. 

Three other points deserve mention in connection with the issue of causal- 
ity. The first is that inflationary particle production would not occur at all 
but for the causal horizon of an inflating universe. Far from being some- 
thing we ignore, causality is at the heart of the effect. The second point is 



that the explicit perturbative computations [H, 11] which support this pic- 



ture were made using the Schwinger-Keldysh formalism |15|. This method is 



manifestly causal: interaction vertices from outside the past lightcone of the 
observation point make no contribution. The final point is that the effects 
of causality are quite evident in the factors of Ht which appear in the result 
(H). They are present, in essence, because two interaction vertices are being 
integrated over the volume of the past light cone (see Fig. 2.), each factor of 
which grows like Ht. 

Those who argue against the possibility of a secular response concede that 
particles continue to attract one another after having fallen out of causal con- 
tact. But they maintain that the lines of force from this effect must be rapidly 
be diluted owing to the inflationary expansion of spacetime. Therefore the 
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Figure 1: The worldhnes of two freely falling particles A and B in co-moving 
coordinates. The interaction A feels from B derives not from some arbitrary 
surface of simultaneity but rather from the portion of its past light cone {B' 
and below) in which B is visible. 



gravitational fields contributed from particles created early during inflation 
should become weaker and weaker. Were this view correct the cumulative 
gravitational potential would not be the integral of a constant — as we have 
argued — but rather some power of the ratio of past to present scale factors, 

rrft'e^(*')-W~^, (5) 
^0 b{t) 

Note that we are assuming a homogeneous, isotropic and spatially fiat back- 
ground geometry, 

ds'^ = -dt^ + e^^^^Ux ■ dx = n^{r]) [drf + dx ■ dx] . (6) 

We have also used the slow roll approximation to evaluate the integral.[| 

There is no question that the argument given above applies for a con- 
formally invariant force field such as electromagnetism. But gravity is not 
conformally invariant, and the gravitational response from early perturba- 
tions approaches a constant rather than redshifting to zero |0. Because it 

iThis means > |6(^)| and also that e''^*) > 1 for iJt > 1. 
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illustrates the crucial distinction between gravitation and electromagnetism 
we will make the argument in detail. Consider an electromagnetic analog of 
the particle creation process in which a homogeneous and temporally con- 
stant current density J\t,x) = j6\ is produced at each point in an inflating 
spacetime. In flat space this would engender a linearly growing electric field 
of the form, E\t,x) = —jtS^, whose interpretation would be the superposi- 
tion over time of a constant electric displacement. 

Things work a little differently during inflation. The curved space Max- 
well equations are, 

(v/^F^^') = J^v^ . (7) 



For our problem only the /i = 3 component is nontrivial. 



dt 

If we assume zero initial electric field the solution is, 

E%t) ^ F'^it) = -J dt' (e^(*')-^W)^ ~ ^ . (9) 

The continuous current produces an electric field which approaches a constant 
during inflation (when h is nearly constant). There is no appreciable buildup 
from previous times because the electric field lines redshift like e~^^.0 

To see what happens without conformal invariance let us first translate 
the electrodynamic problem to the context of a massless, conformally coupled 
scalar. The Lagrangian is, 

C-Mcc = 'Id^^d^ijg^''^ - ^ij'RV^ - J^V^ . (10) 

Specializing its field equation to a constant source J{x) = j in a homogeneous 
and isotropic geometry gives, 

gg^'d.^) -lm = J-^ i e'4eVl = J • (n) 



6 dt \ dt 




^This exercise incidentally illustrates the way in which sources can produce effects even 
when they are not within one Hubble length on some arbitrary surface of simultaneity. To 
see this simply turn the current off at a certain time and watch the subsequent evolution 
of the electric field. Although it decays exponentially it is still nonzero, despite there being 
no current anywhere on the same surface of simultaneity. 
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Assuming no initial field the solution is, 



= -je-^W /* dt'e-'^''^ f dt"e''^'"^ ~ -J- . (12) 
"'0 Jo 2b'^{t) 

Just as in the analogous electrodynamic problem, the conformally coupled 
scalar approaches a constant during inflation because previous contributions 
to it are redshifted. 

We can understand what happens when conformal invariance is absent 
by comparing the response to the same source from a massless, minimally 
coupled scalar. Its Lagrangian is, 

'Cmmc = ~^^(f)^u(pg'"'^/^ - J(t)^f^ ■ (13) 
The same specializations as before reduce the equation of motion to, 

-^d, {V^g^'^d^^p) =J^ -e-''^ {e''<p) = j . (14) 



The solution can again be expressed as a double integral, but now there is 
no outer redshift factor, 

0(t) = -J f dt'e-3^(*') /' rft"e^^(*") ~ -3 f dt'J- . (15) 
Jo Jo Jo hit') 

During inflation h is nearly constant, so ^(t) ~ —jt/h. Hence we learn that 
lifting conformal invariance permits early sources to contribute on an equal 
footing with late ones. 

The scalar comparison we have just given was far from specious. It has 
long been known that, dynamical gravitons in a homogeneous and isotropic 
background geometry obey the same linearized equation as massless, mini- 
mally coupled scalars [Q. It is also well known that a single perturbation of 
the metric does not redshift away during inflation. Instead the (synchronous 



gauge) perturbed metric approaches the form |T^, |T7 



g^^it, x)dx^dx'' — > -dt^ + e^^^^'>-f,j{x)dx'dx^ . (16) 

Boucher and Gibbons prove their "cosmic baldness" theorem by then arguing 
that, since the temporally constant, spatial variation of jijix) must fall off 
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above a certain co-moving Fourier mode, any freely falling local observer is 
eventually unable to sense it. This is true for one perturbation, but what 
happens when there is a mechanism, such as the Uncertainty Principle, which 
generates perturbations with higher and higher co-moving wave numbers? 
There is simply no escaping the fact that the amplitude of 7jj will change 
more and more, and that this change will be visible to local observers. Hence 
the effect is not only real but also secular. 

This brings us to the third objection of principle, that inflationary particle 
production should result in a small fractional increase in the expansion rate. 
The basis of this argument is modeling the stress-energy of particle produc- 
tion as a homogeneous and isotropic classical fluid. If this were correct we 
could compute the gravitational back-reaction using the homogeneous and 
isotropic Einstein equations, 

3b\t) = 3H^ + 8nGp{t), (17) 
-2b{t)-3b\t) = -3H^ + SnGpit) . (18) 

It is straightforward to show that locally de Sitter inflation at Hubble con- 
stant H results in an average energy density and pressure of p = —p = 
H^/IGti'^ per species of massless, minimally coupled scalar.^ Substituting 
b = H and solving the resulting quadratic equation gives the following result 
for the final expansion rate. 




That something must be wrong with this argument becomes apparent 
from the evident fact that it requires the expansion rate to increase. In 
other words, before the steady state solution ( p!9[ ) is attained we must have 
b = —4:7iG{p + p) > 0. But this violates the weak energy condition! Such a 
thing can happen through quantum effects |T^, but the physics is hopelessly 



wrong in this case. A collection of particles, even gravitons, should be a drag 
on the expansion of spacetime, not a super-accelerant to it. 

Closer inspection reveals two flaws in the argument. First, it assumes 
that the distribution of created particles is homogeneous and isotropic. This 
can hardly be so when the density of these quanta is about one per Hubble 

■^For a detailed derivation of this old and well-known result, see Section 4 of flSl. 
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Figure 2: Two loop contributions to the background geometry from gravitons. 
The lefthand diagram is a neghgible constant because the gravitons do not 
interact in a correlated fashion. The righthand one slows expansion by a 
fractional amount which grows like (GA)^(i^t)^ because the gravitons are 
correlated when they interact. 

volume. Gravitons move at the speed of light, and there must necessarily 
be a direction associated with this motion. The one infrared graviton in any 
single Hubble volume cannot give rise to a stress-energy tensor which is even 
approximately isotropic. Pretending otherwise risks the same sort of mistake 
as using the zero average charge density to describe the dielectric properties 
of a ponderable medium which actually consists of an enormous number of 
positive and negative charges. 

The second flaw is ignoring quantum correlations. Fig. 2 presents two 
diagrams that contribute to the quantum gravitational back-reaction on an 
inflating universe. Each has two loops and contributes at order (GA)^. Yet 
the diagram on the left gives a constant that can be absorbed into renor- 
malizing the cosmological constant^ whereas the one on the right contributes 
to the secular slowing of the Hubble constant in expression @ The key 
distinction between the two graphs is that the virtual gravitons on the right 
are correlated when they interact whereas the ones on the left are not. 

It is easy to understand why correlation can make so much difference. 
Consider the traveling wave pulses depicted in Fig. 3. If the stretched string 
upon which they move is exactly linear then the packets pass through one 

^In fact it must be so absorbed in order for the initial condition of inflation with Hubble 
constant H to be satisfied. 
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Figure 3: A simple illustration of how correlation can change an interaction. 

The wave packets on the top have the same sign, so were there a fourth order 
coupling they would repel one another. The wave packets on the bottom have 
opposite signs. With the same coupling they attract one another. 



another with no effect, however we can imagine that there is a nonlinear 
interaction whose contribution to the potential energy of the string goes 
like the fourth power of the displacement. Then it is clear that pulses of 
the same sign and amplitude A repel one another because the superposed 
amplitude of 2A gives a potential energy of +1674"^, which is larger than the 
when they are far apart. On the other hand, pulses with the opposite sign 
attract because their amplitude together is zero, which means zero potential 
as compared with the same 2A^ when far apart. Pulses of the opposite sign 
are a reasonable way (with classical string!) of representing particle- anti- 
particle pairs. Ignoring their affinity for one another — which is what using 
the average stress-energy does — is a terrible approximation. It doesn't even 
get the sign right! 

In considering the effects of anisotropies and correlations it is important 
to distinguish between the bare stress-energy of created gravitons — just the 
hu! per particle — and the stress-energy they develop through gravitational 
interaction with one another. The former is highly anisotropic and also 
highly correlated, which is why one doesn't even get the right sign by using 
its average value to compute the gravitational interaction energy. What the 
various elements of such a source must actually do is to attract one another, 
whereas a precisely homogeneous and isotropic source causes spacetime to 
expand. However, the gravitational interaction at any point is the result of 
superposing the gravitational fields of all pairs that have been created in the 
past light cone of that point. There are no significant quantum correlations 
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between different pairs, and their individual anisotropies tend to average out 
over a long period of inflation. So it should be valid to infer the effect on 
inflation by using (^^-[T^), provided it is the energy density and pressure 
of gravitational interaction that are substituted, and provided that these 
quantities are not themselves computed by using the average of the bare 
stress-energy in the homogeneous and isotropic Einstein equations. 

Of course the correct way of inferring the stress-energy of gravitational 
interaction is from the operator equations of quantum general relativity, and 
precisely that was done to obtain (|]) [§. However, we can understand the 
result by simply forcing the average bare energy density of created particles 
to attract itself as we know that the actual, anisotropic and quantum corre- 
lated source does. Much of the essential physics can even be captured using 
Newtonian gravity to estimate the interaction energy, provided the pressure 
is assumed to follow from conservation. 

Consider locally de Sitter inflation on a manifold whose spatial section is 
a 3-torus. If the physical radius of the universe is initially then its value 
at co-moving time t is, 

r{t) ~ if-^e^* . (20) 

As mentioned before, the average of the bare energy density of infiationally 
produced infrared gravitons is, 

PiR ~ . (21) 

This is insignificant compared with the energy density in the cosmological 
constant (~ H'^/G), and piR is in any case positive. However, the gravi- 
tational interaction energy is negative, and it can be enormous if there is 
contact between a large enough fraction of the total mass of infrared gravi- 
tons, 

M{t) ~ r%t)pi^ ~ He^^' . (22) 

For example, if M{t) was all in contact with itself the Newtonian interaction 
energy would be, 

- -^^ GHe , (23) 

Dividing by the 3- volume gives a density of about —GH^e^^^, which rapidly 
becomes enormous. 
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Of course this ignores causality. Most of the infrared gravitons needed to 
maintain piR are produced out of causal contact with one another in different 
Hubble volumes. The ones in gravitational interaction are those produced 
within the same Hubble volume. Since the number of Hubble volumes grows 
like e^^*, the rate at which mass is produced within a single Hubble volume 
is, 

^ ~ i/= . (24) 

dt ^ ' 

Although most of the newly produced gravitons soon leave the Hubble vol- 
ume, their gravitational potentials must remain, just as an outside observer 
continues to feel the gravity of particles that fall into a black hole. The rate 
at which the Newtonian potential accumulates is therefore, 

dt H-^ dt ^ ' 

Hence the Newtonian gravitational interaction energy density is, 

pit) ~ piR$i(t) ~ -GH^Ht ~ -{GH^f -Ht- — . (26) 

G 

Although this model is very crude it does act in the right sense — to 
slow inflation — and it is secular. The dependence upon coupling constants 
is also correct — (GA)^, characteristic of a 2-loop process like the quantum 
gravitational result (0). This can be understood from the fact that one must 
go to one loop order to see particle production, whereas the interactions 
between these particles — which is the source of the effect — requires another 
order in perturbation theory. 

One other thing that this model gets right is the equation of state, if we 
also assume the relativistic form of stress-energy conservation. When the 
number of e-foldings Ht becomes large, the fractional rate of change of the 
gravitational interaction energy is negligible compared with the expansion 
rate, 

\p{t)\ <^ H\p{t)\ . (27) 
It follows from energy conservation, 

p{t) = -m(p{t)+p{t)) , (28) 
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that the induced pressure must be nearly opposite to the energy density. In 
other words, back-reaction induces negative vacuum energy. Since the key 
requirement is slow accumulation in the sense of relation (PTj), the equation 
of state is really a consequence of the fact that gravity is a weak interaction 
on the scale of inflation. 

Before closing the section we should also comment that one does not re- 
quire a complete solution of quantum gravity in order to study an infrared 
process such as this. As long as spurious time dependence is not injected 
through the ultraviolet regularization, the late time back-reaction is domi- 
nated by ultraviolet finite, nonlocal terms whose form is entirely controlled 
by the low energy limiting theory. This theory must be general relativity, 

with the possible addition of some light scalars. Here "light" means massless 
with respect to H = \Ja/3. No other quanta can contribute effectively in 
this regime. 

It is worth commenting that infrared phenomena can always be studied 
using the low energy effective theory. This is why Bloch and Nordsieck 
[ pUj] were able to resolve the infrared problem of QED before the theory's 
renormalizability was suspected. It is also why Weinberg ||21[ was able to 



achieve a similar resolution for A = quantum general relativity. And it is 
why Feinberg and Sucher were able to compute the long range force due 
to neutrino exchange using Fermi theory. More recently Donoghue has 
been working along the same lines for A = quantum gravity. 



3 Problems with coincident propagators 

As discussed in the previous section, the effect we seek to study has two 
essential features. The first is that infrared virtual particles are continually 
being ripped out of the vacuum and pulled apart by the inflationary expan- 
sion of spacetime. The second crucial feature is that these particles attract 
one another through a weak long range force which gradually accumulates 
as more and more particles are created. The particles we believe were actu- 
ally responsible for stopping inflation are gravitons, and the long range force 
through which they did it was gravitation. However, this is not a simple the- 
oretical setting in which to work. It took over a year of labor to obtain the 
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two loop result (0) — even with computer symbolic manipulation programs 
[^]. Before attempting to extend this feat to include invariant observables or 
stochastic samples one naturally wonders whether there is not some simpler 
theory we could study which manifests the same effect. 

As also explained in the previous section, the prerequisites for inflationary 
particle production are masslessness on the scale of inflation and the absence 
of classical conformal invariance. Massless, minimally coupled scalars have 
these properties — and the lowest order back-reaction from self-interacting 



scalars can be worked out on a blackboard in about 15 minutes [|r^. Of 
course it is not natural for scalars to possess attractive self-interactions and 
still remain massless on the scale of inflation. But we do not need a realistic 
model — that is already provided by gravitation. What we seek is rather a 
simple model that can be tuned to show the same physics, however contrived 
and unnatural this tuning may be. 

The behavior of free, massless and minimally coupled scalars in a locally 
de Sitter background has been much studied [^, ^ ^ Among the 
curious properties of these particles are the absence of any normalizable, de 
Sitter invariant states [p4l and the appearance of acausal infrared singularities 



when the Bunch-Davies vacuum is used with infinite 3-surfaces [25, 26]. The 
feature that concerns us here is the assertion that taking the coincidence 
limit of the propagator gives an ultraviolet divergent constant plus a finite 



term which grows linearly with the co-moving time p7[| . 



tt2 

tA{x; x) = UV + —Ht . (30) 

Although this may not be an observable statement for free scalars we shall 
argue in this section that it has three disturbing consequences for our program 
of adding a X(f)^ self-interaction to obtain a paradigm for infrared quantum 
gravity: 

1. The linear growth derives from the ultraviolet, not the infrared; 

2. Scalars develop a linearly growing mass at order A; and, worst, 

3. The stress-energy violates the weak energy condition at order A. 

The purpose of this section is to demonstrate these problems. We will explain 
how to correct them in the next section. 
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We begin by giving the Lagrangian. Without ordering corrections it is, 
^ = -\d,ct^d,ct^g^^^g - ^A0^v^ + A£ . (31) 
The various counterterms reside in A£, 

-^-6Zd,<pd,<Pg^''V^ - ^6X<P'V^ . (32) 

The ones on the first hne are of order A and will figure in the considerations 
of this section. The ones on the second hne are of order and will not 
concern us further here. 

We are not quantizing gravity. The metric is a non-dynamical background 
which we take to be locally de Sitter in conformal coordinates, 

g^^if], x) = n'ir])r]^, , 1^(77) = --i- = e^* . (33) 

(Recall that A = 3H^.) To regulate the infrared problem on the initial value 
surface we work on the manifold x R, with the spatial coordinates in the 
finite range, —H~^/2 < < H~^/2. We release the state in Bunch-Davies 
vacuum at t = 0, corresponding to conformal time rj = —H~^. Note that the 
infinite future corresponds to rj 0^ , so the possible variation of conformal 
coordinates in either space or time is at most Ax = Arj = H^^. 

Because the spatial manifold is compact, wave numbers have the form, 
k = 27rHn, where n is a vector of integers. Excepting for some completely 
irrelevant zero modes, the free field mode sum has the form 



(pi{r],x) = H^Y.l'^iv, k)e"'-'^a(k) + M*(r/, k)e-'^'^a\k)] , (34) 

k 

where the properly normalized Bunch-Davies mode functions are, 

u{r,, k) = ^ (n-^ + '4-) e-*''' . (35) 



2k \ k J 

If we define At] = rj — rj' the product of the mode function and its conjugate 
can be reduced to the following useful form, 

-ikA-q tt2 
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This is of course the combination which occurs in the propagator. 

We introduce a convergence factor of e""^^ to promote the mode sum for 
the free propagator from a distribution into a well-defined function, 



iA(x; x') 



-ek 



-iklA-ql+ik-Ax 



k 



k 



(37) 



The spatial separation vector and its norm are Ax = x~x' and Ax = ||Aa;||. 
Because the range of conformal coordinates is rather small, it is an excellent 
approximation to represent the mode sum (|37D as an integral. When this is 
done the result is amazingly simple p8| . 



iA(x: x') 



47r2 Ax' - {\Ari\ 



It 



In 



(Ax' - {\Ar]\ - itf)\ + 0(Ar7, Ax) . 



(38) 



Note that the ie term serves as an ultraviolet regulator. It is in fact an 
exponential cutoff on the co-moving momentum. 

We can now demonstrate the first problem by taking the coincidence 
limit. Setting Ar] = Ax = in (pH]) gives. 



iA(x: x) 



47r2 ^2^2 g^2 



In 



2,2 



(39) 



This is the result with a cutoff on the co-moving momentum. For an invariant 
regularization we should really cut off on the physical momentum, fcphys = 
Q~^k. The change is easily made with the replacement e e{HQ)~^, 



iA(x; x) 



H\' r 1 



27tJ 
UV + 



ln(e) + ln(^])l 



2tiJ 



Ht . 



(40) 
(41) 



We are not saying this time dependence is wrong, but its origin is very clear. 
It comes from a logarithmic ultraviolet divergence which is being cut off at 
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+ M 



Figure 4: The scalar self-mass at order A without normal ordering. V denotes 
the 4-point vertex and M stands for the mass counterterm vertex. 

a time dependent point. The terrifically undesirable feature about this is 
that the result can change when one changes the ultraviolet structure of the 
theory, which is precisely the sector of quantum gravity we do not know. 

To see the second problem we compute the one loop scalar self mass- 
squared. The diagrams are depicted in Fig. 4. A trivial application of the 
Feynman rules gives, 

- iM2_i,„p = -'-XtA{x; x) - idm" . (42) 

The ultraviolet divergence must be absorbed with the counterterm, and the 
finite part is fixed by demanding that the scalar be initially massless. There- 
fore the mass counterterm is, 

5m2 = -^UV + 0(A'), (43) 

and the renormalized one loop mass squared is, 

Mi^_i„„p = ±-^H'm. (44) 

This establishes that coincident propagators lead to a linearly increasing mass 
at order A. Of course this is undesirable because gravitons never develop a 
mass. 

To see the third and worst problem, consider the order A (two loop) cor- 
rections to the expectation value of the scalar stress-energy tensor. They 
are given in Figures 5-7. The dominant contribution comes from the dia- 
grams in Fig. 5. Applying the Feynman rules and substituting for the mass 
renormalization (^3|) gives, 

Tl^-' = g,.!^-^[^A{x■x)f-^-^^A{x■x)^, (45) 
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Figure 5: The dominant contributions to the scalar stress-energy tensor at 
order A without normal ordering. V denotes the 4-point vertex and M stands 
for the mass counterterm vertex. 

Note the bizarre nature of the stress-energy tensor with this term. Infla- 
tion actually speeds up. In fact there is a trivial interpretation for what is 
happening: the Uncertainty Principle causes the scalar to wander from its 
classical value of = 0, and that engenders a positive potential energy. It is 
this scalar potential energy that increases the expansion rate. 

Of course the diagrams of Fig. 5 are not the only ones which contribute 
to the stress-energy tensor. Although the others cannot change the [Hty 
terms, they do add important structure to enforce stress-energy conservation 
and cancel the ultraviolet divergences.^ Hence our leading order results for 
the induced energy density (Tqo = —pgoo) and pressure (Tij = pgij), 

p = ^{§^)\iHtr+om}+oi\'), (47) 

P = -^{§^)\iHtr + 0m} + 0iX'), (48) 

""Of particular interest is the conformal counterterm given in the first diagram of Fig. 7. 
In a locally de Sitter background there is no distinction between R and the constant, 12H^. 
However, the distinction does matter in computing the stress-energy tensor and one must 
use the conformal counterterm to cancel overlapping divergences which come from the 
diagrams of Fig. 6. 
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Figure 6: Contributions to the scalar stress-energy tensor at order A which 
involve the derivative vertex D without normal ordering. V denotes the 4- 
point vertex and M stands for the mass counterterm vertex. 

imply that their sum is, 

Note that p + p is not changed if we include the energy density and pressure 
of the bare cosmological constant through the replacements, p — > A/SttG + 
p and p — — A/SvrG + p. This establishes the third and final problem 
with coincident propagators: they result in a violation of the weak energy 
condition at order A. 

4 Covariant normal ordering 

We would like to have a theory in which the ultraviolet does not contaminate 
the infrared, in which massless particles remain massless, and in which matter 
is a drag on expansion rather than a super- accelerant. The purpose of this 
section is to describe a consistent way in which the scalar model can be 
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A 



Figure 7: Contributions to the scalar stress-energy tensor at order A from 
the conformal counterterm C and the counterterm for the bare cosmological 
constant A. 

altered to remove these undesirable features, at least at the lowest orders in 
A. We call the method, covariant normal ordering. For simplicity we shall 
give the result for the original Lagrangian, without its counterterms. Their 
inclusion is straightforward. 

In what follows we assume that the theory has been regulated invariantly 
so that the coincident propagator is a finite scalar functional of the metric. 
We define the normal-ordered product of 0^ as follows. 

Note that one can differentiate either with respect to the field or the coinci- 
dent propagator. 

The trick behind covariant normal-ordering is to implement it at the level of 
the Lagrangian through the replacement: £ — >: jC .. In this way all objects 
derived from the action — such as the scalar stress-energy tensor and the 
scalar field equations — are free of tadpoles. 
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Stress-energy conservation is maintained by taking account of the implicit 
metric dependence of the propagator. One can infer this, formally, from the 
functional integral representation. 



[dmvW) exp [-'- I d'xd,(t>d,(t>g^^^g 



(52) 



Hence the variation which gives the stress-energy tensor produces, 



-2 5iA{y-y') 



-g{x) ^9^"ix 
d d 



X 



{tA{x; x')tA{y; y') + 2zA(x; y)iA{x'- y')},,^, . (53) 



However, we shall modify this scheme in two ways, one necessary and the 
other highly convenient. The convenient modification is that we can drop 
the first of the three terms in ( p3D because it is separately conserved. 

The necessary modification is that one really varies the Schwinger func- 
tional integral I^Uj to obtain the stress-energy tensor. This contains + fields 
which evolve the theory forward and — fields which evolve it back to the 
initial state. Although there is no mixing between these fields, both are min- 
imally coupled to the same metric. Hence the stress-energy tensor receives 
contributions from both terms. This is necessary to make the stress-energy 
tensor real and to make it depend causally upon quantities in the past light- 
cone of the point x^ at which it is evaluated. Two sorts of propagators 
result, 



2A_|_+(a:; x) = 

1 Q-\7])Q-\7]') _ 
47r2 Ax^ - {\Ar]\-iey ~ 8^ 
iA^_{x; x') = 

1 Q~\7])Q-\7]') 



In 



[Ax^ - {\Ar]\ - ief 



47r2 Ax^ - (Ar/ + ze)^ Sn"^ 



In 



(Ax^ - {At] + ie) 



(54) 



(55) 



Of course the ++ propagator is just the same as the Feynman one, iA{x; x'). 
Note also that whereas the free kinetic operator acts upon zA++(x;a;') to 
give a delta function, it annihilates iA+_(x;x'). 
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With the two modifications described above the result is, 




x[dpiA++{x;y)d^iA++{x;y) - dpiA+_{x;y)daiA+^{x;y)] . (56) 



Note that ++ and H — propagators interfere destructively whenever the 
dummy variable of integration, y^, strays outside the past lightcone of the 
observation point x^^. Using the normal-ordered field equations. 



it is easy to see that the stress-energy tensor is conserved. 

Covariant normal-ordering is trivial to use: simply apply the old Feynman 
rules and then ignore any coincident propagators. Including the counterterms 
is straightforward but irrelevant because the first contributions to Sm^, S^, 
6Z and 6X are of order A^. Since the additional terms in the stress-energy 
operator consist of these O(A^) constants times normal-ordered products of 
the fields, the lowest correction to the stress-energy tensor from all except 
the 5 A counterterm are of order A^. This cannot affect the order A^ effect we 
shall compute in the next section. 

Covariant normal-ordering cannot be fundamental because it results in a 
stress-energy tensor that depends nonlocally (but causally) upon the fields. 
However, it does seem to be acceptable if all we seek is a method for tuning 
a simple, scalar model to make it roughly agree with what goes on in the 
vastly more complicated system of quantum general relativity. In particular, 
the model's stability should be enhanced, not endangered, by eliminating the 
tendency of quantum fluctuations to produce super-acceleration at the lowest 
order in perturbation theory. It is also relevant to note that the technique 
reduces to ordinary normal-ordering in the flat space limit. 

5 Back-reaction 

The purpose of this section is to demonstrate that the scalar model deflned 
by covariant normal-ordering shows real back-reaction at the lowest order in 




(57) 
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Figure 8: Contributions to the scalar stress-energy tensor at order with 
covariant normal ordering. V denotes the 4-point vertex and D represents 
the derivative vertex. 



perturbation theory. We begin by evaluating the lowest order contribution 
to the expectation value of the stress-energy tensor. This result is then used 
to compute the expectation value of the invariant expansion operator which 
was defined in a previous paper 0. The conclusion is that back-reaction 
slows inflation by an amount which eventually becomes nonperturbatively 
strong. All orders bounds are given for the strength of the effect. Finally, 
we argue that significant back-reaction would show up as well if stochastic 
samples, rather than expectation values, had been used. 

With covariant normal-ordering the expectation value of the stress-energy 
tensor is much simpler than without. Because there are no coincident prop- 
agators the lowest contribution comes at order from the two diagrams in 
Fig. 8 ITTI. (Of course there is a cosmological counterterm to absorb the 



ultraviolet divergence.) Further, the derivatives on the top vertex of the 
righthand diagram render its contribution subdominant to the lefthand di- 
agram in powers of Ht. So the dominant contribution to the expectation 
value of the stress-energy tensor is —g^u{,x) times. 



: : 



^, / dWW) {[^A++(a;; x')f - [zA+_(x; x')]'} + 0{X') . (58) 
4! Jt'>o ' 

The difference of -|--|- and H — propagators comes from using the Schwinger- 
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Keldysh formalism |]T5|, ^ to compute an expectation value rather than 
an in-out amplitude. This form ensures that the result is real and that it 
depends only upon points x'^ in the past lightcone of the observation point 
x'^. The lower limit of temporal integration at t]' = —H~^ (that is, t' = 0) 
derives from the fact that we release the state in free Bunch-Davies vacuum 
at this instant. 
Although ( 



) was computed in ref. we will go over it in detail. Since 
only the logarithm term of the propagator breaks conformal invariance it is 
perhaps not surprising that the dominant secular effect comes from taking 
this term in each of the four propagators. This contribution is completely 
ultraviolet finite, and its evaluation is straightforward if one goes after only 
the largest number of temporal logarithms. 



4! 



X {in^ [h^ ( 



2103%' 
2732776 




drr 



-In' 



drr Svri In 



H'^ (r^ - (Ar/ + 
H^{Ari'^ - r' 



} 



21n(i7Ar/) +ln(l -x^) 



hi:\HAr]) 



(59) 
(60) 
(61) 



For large Ht the biggest effect comes from the term with the most factors of 
\n{—Hri) = —Ht. That the integrand contributes three such factors follows 
from the expansion. 



ln{HAri) = ln{-Hri') - J2 



n=l 




(62) 



An additional factor comes from performing the integration up against the 
final term in the expansion of the ratio. 



At]^ 

The final result is therefore. 



7/4 



IL 

n'2 



1 



xm' 



2^3' 



{{Htf + 0{Hh^)} + 0{X' 



(63) 



(64) 
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Three points deserve comment before we consider the effect on the in- 
variant expansion observable. First, there is nothing paradoxical about the 
negative sign of the {HtY contribution to expectation value of a positive 
operator. The actual result is dominated by a positive ultraviolet divergent 
constant. It is only after the cosmological counterterm is used to subtract 
this divergence that the ultraviolet finite factor of (Ht)^ dominates the late 
time behavior of the scalar stress-energy tensor at order A^. 

Our second comment is that the negative sign of the {HtY term has a 
simple physical interpretation. As the inflationary expansion rips more and 
more scalars out of the vacuum their attractive self-interaction acts to pull 
them back together. This is the direct analog of the graviton effect we have 
been seeking. 

Our final comment is that the full stress-energy, including the bare cosmo- 
logical constant, obeys the weak energy condition. Our leading order result 
implies. 



Since GA -C 1, the sign of the total energy density is positive, even though 
the scalar contribution is negative. (At least for as long as perturbation 
theory remains valid.) From conservation we see that the sum of the energy 
density and the pressure is positive. 



Pit) +Pit) ^ [liHtf + 0{HH^)] + O(A^) . (67) 



This is the sense in which matter ought to act, so we conclude that covariant 
normal-ordering has succeeded in making scalars behave analogously to the 
vastly more complicated graviton model. 

We have so far been working in a locally de Sitter geometry. To quan- 
tify the back-reaction on inflation we must instead regard de Sitter as the 
background upon which perturbative corrections are superimposed. 




(66) 



(65) 



g;j.u{x) = fl^ 



-1 



(68) 



Hr] 
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Here n = IGttG. The pseudo-graviton field, ipfiuix), is determined, up to a 
diffeomorphism, by solving the Einstein equations, 

Rf,u - ]j^Rg^iu = -^Qfiu + SttGT^^Iq, (p] . (69) 

Note that this need not entail quantizing gravity. It is perfectly consistent to 
suppress dynamical graviton degrees of freedom so that the pseudo-graviton 
field is only an operator through its dependence upon 0, and this is what we 
shall do. One consequence is that Kipniy receives its first nonzero contributions 
at order G. If these result in a secular reduction of the expansion rate — as 
we will see that they do — then we can establish that back-reaction is real 
without needing to consider either higher G corrections to Kipfj,i, or the effect 
of more than a single power of the pseudo-graviton field in the expansion 
operator. 

For the purpose of quantifying back-reaction in this system it suffices 
to use the simplest of the invariant observables previously constructed for 
that purpose 0. We first define a scalar measure of the expansion rate and 
then evaluate it on a geometrically fixed observation point. Our scalar is the 
inverse conformal d'Alembertian acting upon a unit source, 

A[g]ix) ^ (70) 
= Ao{x) + kAi{x) + 0{K^ij^) . (71) 

The zeroth order term can be evaluated exactly, although we shall make the 
slow roll approximation (denoted by an arrow), 

A(x) = (72) 

= _e^^* Tdt'e-^*' r'rff'e^^*" , (73) 
Jo Jo 

Because this is constant, no perturbative shift in the observation point can 
affect the first order result! Therefore, it does not even matter how the obser- 
vation point is geometrically determined although, for the sake of complete- 
ness, we employ zero shift with the following clock function to fix surfaces of 
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simultaneity, 



Ht + 0{KtP) . 



(76) 



The single graviton correction to ^[(y'](x) is simple to evaluate in the slow 
roll approximation H], 



kAi(x) 



K 1 



Q2 



(77) 
(78) 



This particular combination of the pseudo-graviton field happens to be fixed 
by the Einstein equation (|69D, 



kAi (x) 



ttG 



g^^"{x)T,,{x) . 



(79) 



Since we need not consider corrections of higher order in G, the metric in 
this last expression can be set to the non-dynamical background, g^y{x) — > 



^ iv)Vtii^- From the form of the stress-energy tensor (^q) , and our previ- 



ous result for the leading contribution to its expectation value (|6^, we see 
that the slow roll approximation for the expectation value of the expansion 
operator is. 



{n\A[g]ix)\n) 



2H^ 



1 + 



{Hty + 0{H-Y) +0{X\G 



(80) 

It follows that expectation values give the following estimate for the back- 
reacted expansion rate. 



HM = H{1 



X^GA r 



{Hty + 0{HH-') +0{X\G 



(81) 



This is precisely the same result that was previously obtained by studying 
the expectation value of the gauge fixed metric |TT|. Back-reaction is for real. 
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Although we will not compute them here, it is straightforward to estimate 
the strength of higher order effects. Consider a diagram with 2N external 
scalar lines. At L loop order the number of 0^ interaction vertices is, 

V = L + N -I . (82) 

Each contributes a factor of A, so the diagram goes like A^"*"^"^. The number 
of internal propagators is, 

P = 2L + N -2 . (83) 



The largest secular effect comes, as it did for (|64D , from the term where 
each propagator contributes its logarithm part. Since we are computing a 
Schwinger diagram, there will be V cancellations between + and — variations, 
which give the ^-function imaginary part of the logarithm. However, there 
are also V temporal integrations, each one of which can potentially result in 
an extra factor of ln(— if?]) = —Ht. Hence the strongest possible effect for 
the 2iV-point vertex at L loop order is, 

y4~A^+^-i(m)2^+^-2. (84) 

The stress-energy tensor corresponds to = so the dominant contri- 
bution at L loop order is, 

T^^ g,,H^ {XiHtyY'' . (85) 

It follows that perturbation theory breaks down at Ht ~ l/VX. Since A 
is assumed small we see that back-reaction can be studied reliably for an 
enormous number of e-foldings. Note that all the higher point diagrams 
remain perturbatively weak during this entire period, 

lim ^4 ~ A^/2 _ /gg) 
It should therefore be valid to use perturbation theory almost up to Ht = 

i/Vx. 

Finally, we consider the effect of taking stochastic samples of ./l[5](a;) 
rather than computing its expectation value. A procedure for implementing 
this perturbatively was worked out in ref. |T^. What one does is to solve 



29 



the scalar field equations (0) for 0(x) in terms of its initial value and that of 
its first derivative, organized as creation and annihilation operators on free 
Bunch-Davies vacuum, 



A 

(X 



i{x) + - f rfVact(x; x') : (j)^j{x') : +0{X') . (87) 

O Jt'>0 



The free field 4>i{x) was given in expression (|34D , and the retarded Green's 
function is, 



The order G result for the expansion observable comes from substituting 
this solution into (^) and then assigning random (D-number values to those 
creation and annihilation operators in (f)i{x) which have experienced horizon 
crossing by the observation time. 

Since the dominant contribution to the stress-energy tensor is from the 
quartic coupling we can make the replacement, 

\C 

^A,{x) — : <p\x) : . (89) 

The order A contribution to this operator is obtained by replacing all the fields 
(f){x) by the free field (f)i{x). Although this term has zero expectation value 
on account of covariant normal-ordering, a stochastic sample will generally 
be nonzero. We can compute its variance by taking the expectation value of 
its square, 

n) = 2A[iA{x;x)]\ (90) 



b%x) 



Since iA{x;x) grows like H"^ /in'^Ht, we see that the order A contributions 
to kAi(x) fluctuate about zero with standard deviation. 



2H^ 2^3--TT 



{mf + 0{Ht)+0{X\G'')\. (91) 



32.5^3 



For some samples the fluctuation reduces the expansion rate from its classical 
value, for other samples the expansion rate is increased. 

It is important to realize that this order A effect is the same as the one 
previously studied on the nonperturbative level by Linde and collaborators 
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Although its variance does have the stated temporal dependence, what 
an actual observer sees is the effect of the field executing a drunkard's walk. 
Hence the time dependence of the order A contribution to a stochastic sample 
of is not simple. 

The order effect we saw from the expectation value derives from the 
term where three of the fields in : 0^ : are the free field 0/ and the other field 
is the order A correction in (|87|) . The nonlocal character of this term tends 
to wash out its variance [|^, so stochastic samples of it are clustered tightly 
around the mean value (|80|). We conclude that stochastic samples of the 
expansion operator consist of the same secular slowing term we found in its 
expectation value, superimposed upon a genuinely stochastic, random walk 
at order A. Since there is only a small probability for the drunkard's walk 
to exhibit monotonic time dependence, it is possible to distinguish the two 
effects, even if a fluctuation happens to make the order A contribution larger. 
Therefore, back-reaction is still real in the presence of stochastic effects. 



6 Discussion 

In this paper we have employed a simple scalar model to demonstrate that 
there can be significant back-reaction on inflation, even when the effect is 
quantified using an invariant operator measure of the expansion rate and 
even when stochastic effects are included. The expectation value of the in- 
variant expansion operator gives precisely the same result that was previously 
inferred by computing the expectation value of the gauge-fixed metric . 
The situation is more complicated when stochastic effects are included. The 
result in this case is that almost the same secular slowing is superimposed 
upon the perturbative analog of the stochastic effect previously studied by 
Linde and collaborators |^. In both cases back- reaction slows inflation by 
an amount which eventually becomes nonperturbatively strong. 

The scalar model is somewhat contrived in two ways. First, the scalar 's 



natural mass is not zero but rather the scale of inflation, ^jHMpi. Second, the 
covariant normal-ordering prescription of Section 4 results in a subdominant 
contribution to the stress-energy tensor (EBI) which depends causally but 



nonlocally upon the scalar and metric fields. This term plays no role at the 
order we worked, but it is necessary to enforce conservation at higher orders. 
Neither of these features should cast doubt upon the reality of back-reaction 



31 



in quantum gravity. In fact they were imposed upon the scalar model to 
make it more nearly resemble gravity. 

In Section 2 we have also tried to answer the three objections of princi- 
ple sometimes made to a significant back- reaction: causality, redshift, and 
averaging the source. To briefly recapitulate, the gravitational attraction be- 
tween superhorizon particle pairs derives from virtual gravitons which were 
emitted in the past, before each particle exited its partner's causal horizon. 
The same mechanism is responsible for the persistence of gravitational fields 
due to massive objects which have fallen inside the event horizon of a black 
hole. Although electromagnetic fields are redshifted by inflation, gravita- 
tional potentials are not. This derives from the fact that electromagnetism is 
conformally invariant whereas gravity is not. On a concrete level one can see 
it from the theta function term in the retarded Green's function (|88[) which 
is common to minimally coupled scalars and dynamical gravitons. Finally, 
it is not valid to compute back-reaction from the average stress-energy of 
produced particles because this suppresses the quantum correlation between 
produced pairs and because the actual distribution of particles is not uni- 
form. Clumps of energy density attract one another gravitationally whereas 
the perfectly uniform distribution which results from taking the spatial av- 
erage simply increases the overall expansion rate. 

A spinoff of our work is the order A violation ( ^Q]) of the weak energy 
condition when the scalar model is not covariantly normal-ordered. This 
seems to be an analog, on cosmological scales, of quantum violations of the 
energy conditions which have been previously studied on much smaller scales 
0. With observat ions on the current state of the universe not disfavoring 
an equation of state with w < —1 it is worth taking note of models that 
can achieve this. 
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